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ABSTRACT 

Passive microwave satellite observations have frequently been used to observe changes in sea 
ice cover and concentration. Comiso et al. (1989) showed that there may also be a direct 
relationship between the thickness of snow cover (h s ) on ice and microwave emissivity at 90 GHz. 
Because the in situ experiment of Comiso et al. (1989) was limit ed to a single station, the 
relationship is re-examined in this paper in a more general context and using more extensive in situ 
microwave observations and measurements of h s from the Weddell Sea 1986 and 1989 winter 
cruises. Good relationships were found to exist between h s and the emissivity at 90 GHz-10 GHz 
and the emissivity at 90 GHz-18.7 GHz when the standard deviation of h s was <50% of the mean 
and when h s was <0.25 m. The reliance of these relationships on h s is most likely caused by the 
limited penetration through the snow of radiation at 90 Ghz. When the algorithm was applied to 
Special Sensor Microwave/Imager (SSM/I) satellite data from the Weddell Sea, the resulting mean 
h s agreed wi thin 5% of the mean calculated from >1400 in situ observations. 




INTRODUCTION 

The emissivity of sea ice and its associated snow cover is much greater at some microwave 
frequencies than those of seawater. Consequently, passive microwave sensors have proven to be 
very useful for mapping the extent of the sea ice cover. To date, several sensors have been 
launched aboard satellites and have provided a nearly continuous record of sea ice distribution 
since 1972. These include the single band Electrically Scanning Microwave Radiometer (ESMR, 
Nimbus-5 satellite, 1972-1977), the Scanning Multichannel Microwave Radiometer (SMMR, 
Nimbus-7 satellite, 1978-1987) with bands at 5 frequencies (6, 10, 18, 21, and 37 GHz) at dual, 
horizontal (H) and vertical (V) polarization, and the Special Sensor Microwave/Imager (SSM/I, 
Defense Meteorological Satellite Program, 1987-present), a radiometer with frequencies at 19 GHz 

(H and V polarization), 37 GHz (H and V polarization), 85 GHz (H and V polarization), and 22 
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GHz (V polarization only). The spatial resolution of these sensors is approximately 30 km . 

In order to provide in situ measurements which can be compared with satellite observations, 
Comiso et al. (1989) used surface-based sensors to measure local variability in the microwave 
signature of sea ice in the Weddell Sea. Comiso et al. (1989) also examined the effect of snow 
thickness on the emission of microwaves. To do so, they mounted a radiometer ~1 .5 m. above the 
sea ice in an area of substantial snow cover. Radiometric measurements were made of the snow 
surface when it was undisturbed and as the snow was carefully scraped off in layers of 0.05 m. 
Significant increases in brightness temperature (Tb) at 90 GHz were observed when the snow was 
removed, while at other frequencies the signal remained relatively constant (Fig. 1). Comiso et al. 
(1989) demonstrated a nearly linear decrease in Tb with increasing h s at 90 GHz (V) when the 
instrument was oriented 50° from nadir, the approximate viewing angle of the SSM/I. This 
suggested that it might be feasible to estimate the thickness of snow cover within the ice pack using 
satellite mounted sensors. 

The decrease in Tb at 90 GHz with increasing h s has been postulated as caused by greater 
attenuation of radiation emitted from the snow/ice interface by the thicker snow pack. It was argued 
by Comiso et al. (1989) that the absence of the decrease in Tb at lower frequencies was because of 
the smaller influence of volume scattering on Tb, since the wavelengths are large compared to the 
size of the snow particles. Modeling results of the effect of snow cover on sea ice show similar 
behavior. It must be noted, however, that the microwave signal is not only influenced by h s but 
by other characteristics of the snow pack as well, such as the wetness, granularity, density, and 
presence of informalities like ice lenses, as well as the underlying sea ice itself. 
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Unfortunately, the preliminary findings reported by Comiso et al. (1989) were limited to only a 
few measurements made at a single location. The objective of the present study is to re-examine the 
relationship between h s and its associated microwave signature using surface-based data, and if 
possible, to develop an algorithm appropriate for estimating hs from space. The surface-based data 
used in this study were collected during the cruises ANT V/2 (1986) and ANT VDI/2 (1989) to the 
Weddell Sea aboard the R/V Polarstem. Also, satellite TBs are correlated with in situ snow 
measurements where no ship observations were available. 

ANT V/2 (19861 

Observations of both h s and microwave Tg were made during the 1986 Winter Weddell Sea 
Project at 32 out of a total of 37 stations. A ship-mounted radiometer was positioned approxi- 
mately 17 m above the sea ice and measured Tg at frequencies of 6, 10, 18.7, 37, and 90 GHz at 
various angles from nadir. The change in polarization from vertical to horizontal were made by 
rotating the antenna 90°. Only data measured at an angle 50° from nadir were used in the present 
study, because this is the approximate angle at which SSM/I data are collected (i.e., 53°). 

The mean value for h s at each of the 32 stations was determined from about 26 individual 
measurements. The mean value for h s for all stations ranged from 0 to 0.4 m, with 3 stations 
having h s = 0 m, 27 stations with h s <0.25 m and 2 stations with h s >0.25 m (0.26 and 0.40 m). 
The standard deviation of h s ranged from 20-130% of the mean. 

Emissivity (e) at the sea ice surface, calculated as 


TB - Tsky 
8 Ti - Tsky 

(Comiso et al 1989), was highly correlated with Tb in the 1986 data (Fig. 2), where Ti is die 
snow-ice interface temperature and Tsky is the radiance caused by the atmosphere and free space. 
For example, the regression of Tb10 GHz (V and H polarization) against elO GHz (V and H 
polarization) at the same station yielded a correlation coefficient of >0.99 (Fig. 2A). This was also 
true at 90 GHz (V and H polarization) (Fig. 2B). 

A plot of e90 GHz (V) versus s8.7 GHz (V) (Fig. 3) illustrates that the latter is reduced at h s =0 
m (i.e., the three data points in the lower right quadrant). When snow was present on the surface 
of the sea ice, the signal at 18.7 GHz remained approximately constant, while the micro-wave 
signal at 90 GHz varied from 0.75 to 0.98 (Fig. 3). This relationship was shown previously by 
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Comiso et al. (1989) who attributed the variability at 90 GHz to differences in either h s or 
properties of the snow cover. A similar relationship also was apparent when the 90 Ghz (s and 
Tb, V and H polarization) band was plotted against 10 or 37 GHz bands, but not the 6 GHz band 
which exhibited more scatter. 

Plots of h s versus e or Tb at 6, 10 (Fig. 4A), 18.7, and 37 GHz remained relatively flat (no 
slope), indicating that there was little sensitivity in the passive microwave signal at these frequen- 
cies to variation in h s within the snow pack. The obvious exception was the case where h s =0 m 
which exhibited much lower values for both e and Tb On the other hand, the 90 GHz (V) fre- 
quency exhibited much greater variability with respect to h s , although there was no significant 
correlation between the signal at 90 GHz and h s (Fig. 4B). 

However, when stations where the standard deviation of h s was greater than 50% of the mean 
were removed from the analyses, the relationship between e90 GHz (V) and h s improved dramati- 
cally (Fig. 5). Two types of outliers in Fig. 5 can be identified; stations where h s = 0 m (n = 2) 
and stations where h s > 0.25 m (n = 2). When those stations were omitted, a good linear relation- 
ship between h s and the 90 GHz (V) frequency was obtained (Fig. 7). The equations and 
goodness of fit of the least-squares linear regression between h s and s90 GHz (V) and Tb 90 GHz 
(V) are: 


e90 GHz (V)= 1.08 -0.019 *h s (cm) (n = 10, r 2 = 0.78) (1) 

Tb 90 GHz (V) = 290.24 - 5.21 * h s (cm) (n = 10, r 2 = 0.77) (2) 

Measurements with vertical polarization were more strongly correlated with h s than those with 
horizontal polarization, which had higher (by -10%) standard deviations. The number of stations 
where measurements at vertical polarization were taken was higher as well. The linear regression 
of the h s data from Comiso et al. (1989) was 

Tfi90 GHz (V) = 257.89- 1.06 * h s (cm) (n = 4, r 2 = 0.99) (3) 

and is shown in Fig. 6 along with die data used to generate Eq. (2). Equation (2) has a much 
steeper slope than Eq. (3), suggesting that predictions of h s made using Eq. (2) will be more 
sensitive to changes in T90 GHz (V). Also, if the measured Tb is below 250°K, Eq. (2) will 
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predict lower values of h s , while at Tb above 250°K Eq. (2) will give higher values. At T 250°K, 
both algorithms produce a similar value for h s (-0.08 m). 

If the relationship between h s and microwave signal is expanded to include two frequencies, 

90 GHz (V) and either 18.7 GHz (V) or 10 GHz (V), the regression improves further still 


e90 GHz (V) - elO GHz (V) = 0.20 - 0.0227 * h s (cm) (n=l 1 , r 2 =0.92) (4) 

Tb90 GHz (V) - TBlO GHz (V) = 55.56- 6.20 *h s (cm) (n=12, r 2 =0.92) (5) 

e90 GHz (V) -818.7 GHz (V) = 0.13- 0.0197 *h s (cm) (n=12, r 2 =0.89) (6) 

TB90 GHz (V) - Tb18.7 GHz (V) = 35.91 - 5.43 * h s (cm) (n=12, r 2 =0.90) (7) 


The advantage of this approach was that, in contrast to Eqs. (1) and (2) where snow-free stations 
were omitted, stations where h s = 0 m fitted the resulting regression well (Fig. 7). Thus, the 
algori thm s shown in Eqs. (4-7) apply to all cases where h s < 0.25 m. This restriction is because 
of the fact that at higher values of h s , the penetration depth of radiation at 90 GHz is smaller than 
the snow depth itself (Hall & Martinec 1985). 

Only Eqs. (6) and (7) can be applied to SSM/I satellite data to generate estimates of h s because 
unlike SMMR, SSM/I does not have a 10 GHz channel. SSM/I does not have a 90 GHz channel 
either, but it is expected that the differences between the data generated from the 85 GHz frequency 
of the SSM/I and the 90 GHz field data will be small enough to allow substitution of the 85 GHz 
SSM/I data into Eqs. (6) and (7). 

Foster and Chang (1993) developed a simple algorithm to predict h s on land using the 
microwave brightness temperatures at 18 and 37 GHz. Their equation was 

h s = 1 .59 * [Tb 1 8 GHz (H) - Tb37 GHz (H)] (8) 

However, this algorithm provided a poor fit to the sea ice data, perhaps because of the different 
emissivities of land and ice. 

ANT VTII/2 U9891 

Du ring the Weddell Sea Experiment, h s and surface based radiometer measurements were taken 
simultaneously at 18 stations at an angle 50° from nadir from September 16 to October 6, 1989. 
Unlike 1986, only 10 GHz (V and H polarization) and 85 GHz (V and H polarization) frequencies 
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were measured. Although more detailed information of the snow cover was reported (grain size, 
vertical structure of the snow layer, presence of ice lenses, etc.), mean h s at each station was 
determined from, at most, four measurements within the radiometer footprint. Snow depths 
ranged from 0.08 to 0.61 m and were generally less variable than in 1986. Most of the thick snow 
cover was observed in the western Weddell Sea where the ice cover is generally thicker. Standard 
deviations for a specific station could be >100% of the mean hg, but were generally <40%. 

Plotting the 85 GHz versus the 10 GHz frequency conformed to expectations; Tb at 10 GHz 
was fairly constant, while the signal at 85 GHz was highly variable. Also, Tg and e for the 10 
GHz channel were highly linearly correlated (r2>0.99) (Fig. 8A). However, Tb and 8 were 
poorly correlated for the 85 GHz band for both H and V polarization (Fig. 8B). This observation 
was in direct contrast to measurements made in 1986 which showed a highly linear relationship 
(Fig. 2) and suggests that the sky calibration of the radiometer in 1989 was not as good as those in 
1986. However, some of the 1989 data conform to the relationship between Tb and e found in the 
1986 data. When the line of best fit from 1986 (Fig. 2) is plotted with the 1989 data, some of the 
points lie along the regression (Fig. 8B), suggesting that the calibration for these points may be 
correct. It is interesting to note that all of the suspect data lie to the right of the 1986 regression 
line. 

Because of the calibration problem, h s was poorly correlated with the quantity (e85-el0) GHz 
(V) in the 1989 data set (Fig. 9). However, if the tine of best fit of this relationship from the 1986 
data is plotted with the data from 1989, it can be seen that some of the 1989 data points lie on or 
near the 1986 regression tine. Interestingly, the 1989 data points that conform to the 1986 relation- 
ship between Tb85 GHz and s85 GHz also conform to the 1986 relationship between h s and the 
quantity (e85-10) GHz (V) (see circled solid black data points on Figs. 8 and 9). There were four 
other stations in 1989 which conform to the 1986 relationship between Tg 85 GHz and e85 GHz 
but do not fit the 1986 relationship between hs and the quantity (e85-10) GHz (V). Therefore, the 
fact that Tb and 8 at a particular station relate as expected cannot be used as a proxy for “good 
data” for the snow depth algorithm. On the other hand, the disagreement of those data with the 
algorithm could be caused by experimental or weather conditions as well. For example, the header 
of these data files contained comments such as: “probably moist snow surface,” “radiometers may 
not be warm enough yet,” “cloud broken/sun in and out station,” and “near white out start of 
angular scan, slow motion work to be done at the end”; conditions which would very likely affect 
the microwave signature. 
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Comparisons with SSM/I data 

Eqs. (6) and (7) were used to predict hs in the Weddell Sea from SSM/I imagery. The mean 
value for hs estimated from SSM/I data was nearly identical to that calculated from >1400 in situ 
measurements. However, the frequency histograms exhibited different distributions, with the in 
situ data exhibiting a lognormal distribution and the estimates from SSM/I tending toward a more 
normal distribution (Fig. 10). The differences in the two distributions is not surprising, however. 
The Central Limit Theorem states that the sampling distribution of the mean approaches the normal 
distribution as the sample size increases, regardless of the shape of the distribution of the raw data. 
In our case, the SSM/I averages the snow depth over an area of 625 km 2 , which is essentially the 
same as increasing the sample size, resulting in a normal distribution. Despite the differences in 
distribution, the agreement between the means of the two data sets suggests that the algorithm is 
able to provide a reasonable estimate of h s from SSM/I satellite data over a large area. 
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FIGURES 







Snow thickness (cm) 


Figure 1. Brightness temperature at 90 and 37 GHz as a function of snow depth. Data was 
collected using a microwave radiometer mounted ~1.5 m. above the sea ice in an area of 
substantial snow cover. Radiometric measurements were made as the snow was carefully scraped 
off in layers of 0.05 m. 
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690 GHz 

Figure 2. Brightness temperature versus emissivity at frequencies of (A) 10 GHz and (B) 90 GHz. 
Data was collected in 1986 as described in Fig. 1 with the exception that the radiometer was 
mounted on the RV Polarstem ~17 m above the ice surface. 


12 




0 10 20 30 40 50 



10 20 30 40 50 

Snow depth (cm) 


Figure 4. Snow depth versus emissivity at frequencies of (A) 10 GHz and (B) 90 GHz. Data was 
collected as described in Fig. 2. 
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Figure 5. Snow depth versus emissivity at 90 GHz after removing data points where the standard 
deviation of the snow depth was <50% of the mean. Circled data points indicate where h s = 0 and 
hs > 0.25 m. Data was collected as described in Fig. 2. 


15 


260 


N 


250 L 


240 L 


Comisoet al. (1989) 


O 230 L 


o 

05 

CD 


220 


210 


200 



5 10 15 

Snow depth (cm) 


20 


Figure 6. Least-squares regression of brightness temperature at 90 GHz versus snow depth (closed 
squares). Data points where the standard deviation of the snow depth was <50% of the mean, and 
where hs = 0 and hs > 0.25 m were also omitted. The best fit line of the data from Comiso et al. 
(1989) is shown for comparison (open diamonds). Data was collected as described in Fig. 2. 
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Figure 7. Snow depth versus the difference in emissivity at 90 GHz and 18.7 GHz. Note that data 
points where h s = 0 m have now been included. Data was collected as described in Fig. 2. 
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Figure 8. Brightness temperature versus emissivity at frequencies of (A) 10 GHz and (B) 90 GHz. 
Data was collected in 1 989 as described in Fig. 1 with the exception that the radiometer was 
mounted on the RV Polarstem ~ 1 7 m above the ice surface. The line in (B) is the least squares fit 
to the 1986 data shown in Fig. 2B. Circled data points in (B) correspond to circled data points in 
Fig. 9 and indicate data which may be salvageable. 
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Figure 9. Snow depth versus the difference in emissivity at 85 GHz and 10 GHz. The line is the 
least squares fit to the 1986 data shown in Fig. 7. Data was collected as described in Fig. 8. 


19 




Snow thickness (m) 


Figure 10. Comparison of the relative frequency distributions of snow depth as determined from in 
situ data and from estimates made using SSM/I data as input to Eq. (6). Note that the means for 
the two methods are almost identical. Differences in the distributions are explained in the text. 
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